Calculating energy shifts in terms of phase shifts 
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To clarify the relation of energy shifts to scattering phase shifts in one-body and many-body 
problems, we examine their relation in a number of different situations. We derive, for a particle in 
a container of arbitrary shape with a short-range scattering center, a general result for the energy 
eigenvalues in terms of the s-wave scattering phase shift and the eigenstates in the absence of the 
scatterer. We show that, while the energy shifts for a spherical container are proportional to the 
phase shift over large ranges, those for a cubic container have a more complicated behavior. We 
connect our result to the description of energy shifts in terms of the scattering T-matrix. The 
general relation is extended to problems of particles in traps with smoothly varying potentials, 
including, e.g., the interaction of a small neutral atom with a Rydberg atom. We then consider the 
many-body problem for particles with a two-body interaction and show that the energy change due 
to the interaction is proportional to an average of a generalized phase shift that includes the effects 
of the medium. Finally, we discuss why, even though individual energy levels are very sensitive to 
boundary conditions, the energy of a many-body system is not. 

PACS numbers: 03.65.Nk, 67.85.-d, 31.15.-p 



I. INTRODUCTION 

The art of calculating energies of many-particle sys- 
tems in terms of scattering phase shifts has a long history, 
dating back to Lenz [![, who derived the low-energy ef- 
fective interaction between particles in terms of the scat- 
tering length for two-body scattering, and to Beth and 
Uhlenbeck , who expressed the second virial coefficent 
for a classical gas in terms of the phase shifts for two-body 
scattering. Such a procedure is convenient, since it sep- 
arates the physics of the short-range scattering process 
from that of single-particle motion in the external poten- 
tial confining the particles. There are many variants of 
the formalism, including the pseudopotential approach of 
Huang and Yang || . 

Further applications, particularly to systems in finite 
geometry, include the early example of calculation of en- 
ergy shifts of electronic states due to impurities in a metal 
[H, and studies of ultracold atomic gases confined in 
traps. Yet further examples are calculations of Rydberg 
molecules the calculation of the energy levels of two 
interacting particles in a harmonic trap [6, tJ, and the 
construction of generalizations of the Gross-Pitaevskii 
equation to include interactions beyond the usual pseu- 
dopotential proportional to the s-wave scattering length 
between the particles [1, Q . From a quite different per- 
spective, Liischer [l(| derived a formula for the energies 
of bound states of two particles in a box (with peri- 
odic boundary conditions) that is enjoying use in parti- 
cle physics to extract scattering properties of elementary 
particles from lattice calculations of energies [ll[ . 

Our primary aim in this paper is to clarify the relation 
between energy shifts arising from interparticle interac- 
tions and the scattering phase shifts. We begin by con- 



sidering the energy levels of a particle in a spherical box 
in the presence of a scatterer at the center of the box. 
For simplicity, we focus on energies low enough that only 
s-wave scattering from the scatterer need be taken into 
account. In the presence of a short-range scatterer, the 
energy shifts are linear in the s-wave phase shift, 8, and 
not a trigonometric function of 8, such as the real part 
of the scattering amplitude. However, in the case of a 
particle in a more generally shaped container, with more 
complicated boundary conditions on the wave functions, 
we find that the energy shifts depend implicitly on tan 5 
(see Eq. (|14p ). In particular, as we shall see explicitly for 
a particle in a cubic container, while the energy shifts are 
approximately linear in 8 for small 8, the slope of the shift 
varies from state to state. In the many body problem, 
as we illustrate with a simplified separable interaction, a 
delta-shell potential, the energy shift is linear in an effec- 
tive scattering phase shift that takes into account effects 
of the medium. 

We begin, in Sec. II, by reviewing the energy shifts 
of a particle in a spherical container. We then go on in 
Sec. Ill to derive a general expression for the energy of a 
single particle moving in a container of arbitrary shape 
in the presence of a scattering center. The energy is ob- 
tained as an implicit function of the s-wave phase shift of 
the scatterer and the eigenfunctions of the one-particle 
problem in the absence of the scatterer. We apply the 
general result to the case of a cubic container in Sec. IV. 
In Sec. V we derive results for the cubic container from 
a T-matrix approach, and in Sec. VI describe applica- 
tions to situations such as atoms in traps and Rydberg 
atoms in which the particle moves in a potential that 
varies smoothly in space, rather than in a container with 
infinitely repulsive walls. Section VII considers applica- 
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tions to the many-body problem, and discusses how the 
formalism may be generalized to take into account the 
effects of the medium on two-body scattering. We take 
H = 1 throughout. 



II. ENERGY SHIFTS IN A SPHERE 

To set the stage, we first review the relation of the en- 
ergy shifts and s-wave scattering phase shifts for the fa- 
miliar problem of a particle of mass M interacting with a 
short range potential v(r) in a hard walled spherical con- 
tainer whose radius, R, is much larger than the range, b, 
of the potential. At low energies, scattering is predom- 
inantly s-wave, even if the potential is not central, and 
may be parametrized in terms of the s-wave phase shift, 
5. Consequently, at low energies, eigenstates of the en- 
ergy are also eigenstates of the angular momentum, and 
we shall characterize states by the value of the angular 
momentum. In this limit, states other than s-wave ones 
are unaffected by the scatterer. For v(r) = 0, the single- 
particle eigenenergies for s-wave states are E v = k^/2M, 
with k v = vtt/R, and v = 1,2,3, ... , with eigenfunctions 



1 sva.(k v r) 



(1) 



y/2irR r 

In the presence of the scatterer, the eigenfunctions at 
points beyond the range of the potential become 



sin(k' v r + 5(k)) 



(2) 



where A\ is a normalization constant. The hard wall 
boundary condition implies that 



or 



klR + 6(K)=VTT, 



cot{k' u R) = - cot S(k' v ). 



(3) 



(4) 



The shift in the wavevector, Ak u = k' u ~k u — —8{k' v )/R, 
is linear in the phase shift. Without loss of generality we 
may take 6 to lie in the range —tt/2 to tt/2, since an in- 
crease of 5 by 7r takes one from a given state in the sphere 
to the next higher state. Since beyond the range of the 
potential ipl(r) is a solution of the Schrodinger equation 
with zero potential, the corresponding eigenenergies are 
E v = k' v 2 /2M , and the shifts in the energies due to the 
interaction are 



AE„ = 



k' v 
2M 



Mk„ 



1 - 



SJK) 

2k v R 



IVv(o)| 2 



(5) 



The term quadratic in S may be neglected if S <C 1, 
since k v R > 1. Even for resonant scattering, 5 ~ 7r/2, 
the quadratic term may be neglected for states with many 
nodes (k„R ^> 1). The energy shift is then given by 



AE V ps - 



(6) 



Thus in these general situations, the energy shift is linear 
in the phase shift (and not, e.g., the real part of the scat- 
tering amplitude ~ e lS sin S). The energy of a state must 
be determined self-consistently, since S depends on en- 
ergy, but if R~ 1 (\d5/dk\) -C 1 one may replace k' u by the 
wave number of the state in the absence of the scatterer, 
k v . 

In the limit k — ¥ 0, away from a scattering resonance, 
S — > —ka s , with a s the s-wave scattering length. Then 
Eq. ^ gives 



AE{k) 



2na i 



1-0(0)1- 



(7) 



III. 



S-WAVE SCATTERING IN AN 
ARBITRARY CONTAINER 



We turn now to establishing, for a particle of mass 
M interacting with a short range central potential in a 
container of arbitrary shape, the relation between the en- 
ergy shifts of the particle and the s-wave scattering phase 
shift. We consider a particle interacting with a scatterer 
located at some point within the container, which we take 
to be r = 0. We shall take the potential to be spheri- 
cally symmetric and short-range, vanishing for r > b, 
and we assume that the characteristic dimension of the 
container, A, is much larger than the range of the po- 
tential, b. We shall consider energies, E, sufficiently low 
that only s-wave scattering is important. We denote the 
eigenfunction ipk{ r ) for energy E = k 2 /2M by ipf{r) for 
r < b and for r > b. For distances r <C A, ^(r) 

may be expanded in terms of spherical waves as 



sin(fcr) 



tan S 



cos(kr) 



kr kr 

oo / 

+ X! a imji(kr)Yi m (6,(j)) 

l—l m— — l 



(8) 



where j; is the spherical Bessel function, the a; m are ex- 
pansion coefficients, and A2 is a normalization constant. 
In the sum, terms proportional to the spherical Neumann 
function n/, which is singular for r — > 0, have been omit- 
ted, since we assume that only s-wave scattering is im- 
portant, and consequently the wave function for other 
partial waves must be finite for r — ¥ 0. Thus asymptoti- 
cally, in the region b < r <C 



tan 5 
kr 



+ l + 0(kr) 



(9) 



To derive an expression for the energy of a state, we 
now find a second expression for the wave function, by 
expressing it in terms of the single particle eigenstates 
{(f> n (r)} for v(r) = 0, with energy E n as follows. The 
single particle Green function for v(r) — describing 
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propagation of a particle from the origin to r is given by 



E n — E 



This satisfies the equation, 



1 



— (V' + k')G fc (r,0) = -*(r), 



(10) 



(11) 



which shows that Gfc(r, 0) satisfies the Schrodinger equa- 
tion for V'fe > ( r ) when r ^ 0. Thus outside the range of 
the potential, Gfc(r, 0) satisfies the Schrodinger equation 
and in addition satisfies the same boundary conditions 
as the wave functions 4> n at the walls of the container for 
any boundary condition that is linear in the wave func- 
tion. Equations (fTTj) implies that Gfe(r, 0) = M/2nr for 
One can in fact show that for r > b 



2irA2 tan 5 
Mk 



G fe (r,0), 



by following Ref. [T^. The function 



Xfe(r) = ^T(r) - 



2irA2 tan <5 
Mk 



Gfc(r.O) 



(12) 



(13) 



is regular everywhere inside the container and is an eigen- 
function with eigenenergy E for v(r) = 0. However, un- 
less E = E n , this is impossible, since the (f> n form a 
complete set, and therefore Xfc( r ) — 0- (F° r discussion of 
the more complicated case when E — E n , see Ref. jloj|.) 

Matching the singular and nonsingular parts in 
Eqs. © and (TT2"j) for kr -C 1, we obtain an expression for 
the energy of the state: 



k 

— cot S(k) = lira 
47T r- >-0 



E 



1 ^(r)C(O) 1 
2M - £ 47rr 



(14) 



since non s-wave contributions vanish in the limit r — > 0. 

Equation (TT4")) is one of the main results of the pa- 
per, and determines the eigenenergies in the presence 
of the potential. When the potential strength tends to 
zero, E approaches an unperturbed energy £ ; accord- 
ingly Eqs. (fTUj) and (TT2"j) implies that V'fe(r) approaches 
the state 



E n =S 



(15) 



where the sum is over all possible degenerate states of 

— 1/2 

energy £, and C = (J2e =£ \^n(0)\ 2 ) i s tne normal- 
ization constant. By construction, </>§(r) has a nonzero 
s-wave component, i.e., </>§(0) 7^ if not all <^>*(0) 
equal zero. The remaining states </>f (r) in the subspace 
spanned by the states with E n = £ vanish at r = 0, 
as we now show. Since these states are energy eigen- 
states in the absence of the scatterer, we may write 
<^t( r ) = Se„=£- c n4>n{^)- We now evaluate the inte- 
gral / dr (<^>§(r))* </>§(r), which by orthogonality is zero 



if s 7^ s. Inserting Eq. (|15p and the expansion for 
</>§(r) in the integral, one finds J dr (</>|(r))* </>§(r) = 
CJ2 En =£ c nM0) = C^|(0) = 0. Therefore the </>f(r) 
for s 7^ s vanish at the origin and these states will be 
unaffected when the scatterer at r = is introduced. 

For k — > 0, away from a scattering resonance, the left 
side of Eq. (fl"4| becomes — l/4-7ra s , and so for a s small 
compared with a characteristic dimension of the con- 
tainer, E must be very close to an unperturbed eigenen- 
ergy E n . Changing the order of taking the limit r —> 
and summation in Eq. (|14|) . we obtain 



AE = E 



(16) 



thus the result ([7]) found in a sphere holds for more gen- 
eral geometry. 

We now summarize the salient results of this section. 
When the energies of states of a particle in the container 
in the absence of the scatterer are degenerate, one may 
construct one state, Eq. (|15l) . that does not vanish at 
r = 0, while all other orthogonal states in the degen- 
erate manifold vanish at r = 0. For weak scattering, 
\S(k)\ <C 7r, the energy shift of the first state is given by 
Eq. (fTBj) while the shift of the other states vanishes. In 
our derivation of these results, we have taken into ac- 
count degeneracy of the energy levels in the absence of 
the scatterer. For stronger scattering, it is necessary to 
solve Eq. (Till) . 

The energy shift AE is reproduced by the expectation 
value of the effective interaction 12], 



Kff(r) 



2na t 



(17) 



in the unperturbed state, 4> n . In terms of the interaction 
of two particles of mass m, the coefficient of the effective 
interaction is 47ra s /m. 

One should note that the effective interaction (IT71) is 
not useful for calculating the perturbed wave function, 
^>(r), outside the range of the potential in terms of a 
boundary condition at the origin. Here one needs rather 
to employ the pseudopotential Q 



27rtan(T d 



in the Schrodinger equation 



2M 



-^7 + V P s(*)H»(r) = ^>(r) 



2M 



(18) 



(19) 



The solution is clearly Eq. ([5} . The present discussion is 
readily generalized to higher angular momentum, with a 
pseudopotential involving higher derivatives of the wave 
function IRa. 



IV. ENERGY SHIFTS IN A CUBIC BOX 

We now apply Eq. (IT4l to the specific example of scat- 
tering by a spherically symmetric potential centered in a 
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cubic box of side L, with periodic boundary conditions. 
This geometry, used in the calculations of Refs. [ToL flllj . 
is more difficult to calculate analytically than in a spher- 
ical container, since the boundary is incompatible with 
the rotational symmetry of the central potential. As one 
sees explicitly by direct construction of the solutions of 
the Schrodinger equation, all partial waves are mixed in 
determining the eigenenergies [loj . 

In the absence of the scatterer, the states are simply 



6p(r) 



1 



(20) 



where p = {p x ,Py,Pz} = {2it/L){l x , l y , l z } with 
l x ,l y ,l z — 0, ± 1, ± 2, . . . . As before, we consider only 
the s-wave contribution to the interaction, and from 
Eq. ((Ill) find [Io| 



— cot 8(k) 

47T 



lim 

r->0 



1 



L 3 p 2 
P 



1 

Airr 



(21) 



A relation similar to Eq. (|21[) is more useful for simu- 
lations on a lattice II ill: 



4tt 



cot S(k) 



1 



|p|<A 



1 



p2 _ fc2 



A 

47r' 



(22) 



where A is the momentum cutoff on the lattice. To derive 
Eq. (|2"2"|) , we note that in the limit r — > the terms in the 
sum in Eq. (|2ip for |p| > A and k <C A can be written as 



1 

L 1 



|p|>A 



d 3 p 



o«pr 



d 3 p 



o«pr 



= — - — . (23) 
4?rr 



i2;T) :; ,fi J n (2tt) 3 p 2 

A 

47T 

For small 8, the energy shift obeys Eq. (|16l) . Whereas 
the shift in k for a sphere is linear in — 8, for a cube it is 
only approximately linear in — S, as we see from Fig. (fTJ), 
which gives a plot of k (panel (a)) and E (panel (b)) vs. 
—8 for the first few eigenstates, calculated from Eq. (f2"2"|) . 
When 8 is small (except in the ground state), 



k[L/2-K « kiL/2-iT- a kl S, 



(24) 



where E = kf /2M is the eigenenergy in the absence of 
the scatterer, and from Eq. (|2"2")l . a kl = Dk t /(kiL) 2 , with 
Z?fc, = J2 p dp 2 ,k 2 the degeneracy of the level E. Therefore 
Ak' l L/2n « -8D kl /{kiL) 2 for small 8, where Ak[ = k[ - 
fe;. The linear relation is valid only for small 8, but with 
a slope which varies with the eigenstate of the cubic box. 

The different dependences of k and E on 8 for cubic and 
spherical boxes are due solely to the different boundary 
conditions. The information on the boundary conditions 
is encoded in Eq. (TT4"]) through the unperturbed eigenval- 
ues and eigenfunctions, while the effect of the short-range 
potential is taken into account via 8; in an arbitrary con- 
tainer, the energy change is a function of cot 8 but the 
exact form of the function is determined by the boundary 
conditions. 




-n/2 -tt/4 



tt/4 n/2 

-8 



E/(2tt 2 /ML 2 ) 




FIG. 1: (a) Wavenumbers k and (b) eigenenergies E vs. —8 
in a cubic box. The curves are calculated from Eq. (|22p with 
A = 200tt/L. 



T-MATRIX APPROACH 



With a view to generalizing results to situations other 
than dilute gases, where binary collisions are the most 
important interaction process, it is useful to derive the 
result (j2"2"j) in terms of the T-matrix. The formalism 
may then be applied to calculate properties of interacting 
many-body systems, when the elementary excitations are 
quite different from free particles, e.g. liquid 3 He at low 
temperatures. Equation (|14p was derived by inspecting 
the Schrodinger equation with boundary conditions at 
both small and large r. The T-matrix formalism we give 
below is equivalent to solving the Schrodinger equation 
in the one-body problem and its analogue in the many- 
body problem to be discussed later provides a convenient 
way to treat the interparticle interactions. 
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The T-matrix is defined symbolically as 



Using 



T(E) =v + v 1 T(E) = v + v ' 



E-Hn 



E- H 



(25) 



where Hq is the free particle Hamiltonian, v is the po- 
tential operator, H = H + v, (E — H)^ 1 is the particle 
propagator, and E the energy; this equation shows that 
the T-matrix has poles in E at the eigenenergies of the 
full Hamiltonian. To derive Eq. ([22]) we rewrite ([25]) as: 



T~ 1 (E) = V 1 



1 



(26) 



E - H ' 

with the on-shell T-matrix, at E = k 2 /2M + in, given by 

Mk 

T-\k 2 /2M + in) = (cot S(k) - i). (27) 

2n 

For simplicity we assume a short ranged contact inter- 
action, U6(r) with a momentum cutoff, A. Then in free 
space Eq. (|26|) becomes, 



d 3 p 



T-\E)=U- 1 -[ p 
Jo 



1 



i) 3 E-p 2 /2M 
Similarly in the cubic box, 



L 3 ^ E~p 2 /2M' 

|p|<A y ' 



(28) 



(29) 



Subtracting Eq. ([29]) from Eq. ([28]) , evaluating both sides 
at E = k 2 /2M + ir), where k 2 /2M is now an eigenenergy 
in the box, using Eq. ([27]) and the vanishing of Xi^ x at 
k 2 /2M, we find 



k I txtiA ^ f A d3 P 1 



-y 

r,3 



L 3 k 2 + in — p 2 

p|<A 



(30) 



The integral on the right in the limit k <C A is — (A + 
ik)/4n. Dropping the in in the sum, we recover Eq. (1221) . 

To see that the general formalism leads to Ak n R = —6 
for a particle in a spherical box (cf. Eq. ([3])), we apply 
the T-matrix formalism. For s-waves in a spherical box, 
the analog of Eq. ([2"9l is 



E 



2irR ^ E~p 2 /2M' 

P<l"'m.x/fl 



(31) 



where the p's take on the values vtt/R, with v = 
1, 2, 3, . . . , z^ max . By the same argument as for the cu- 
bic box, we derive the result 



2 . 2k 2 ^ 
k cot 5 = A + 2z/ max — > 



1 



p<7ri/ max /H 



k 2 — p 2 



(32) 



E 



we obtain 
cot 5(k) = - cot(kR) 



1 

kR 



= cot z, 



2_ ( AR 
kR \~ 



(33) 



(34) 



Thus to recover the correct relation, AkR = —5, we need 
to choose the continuum cutoff A to be n{v ma x + l/2)/i?, 
i.e., half way between successive v values. This subtlety 
shares a similar origin as in the transformation of a sum 
of l/n into an integral leading to Euler's constant 7 « 
0.577. An alternative way of dealing with this problem 
is to introduce a smooth, rather than a sharp, cutoff. 

To further illustrate the T-matrix approach, we con- 
sider the delta shell potential 



v ( r ) = {va/Anrl)5(r - r ), 



(35) 



which has the advantage that T-matrix can be con- 
structed explicitly without the need for a cutoff. In the 
continuum limit, where we use as a basis the s-wave func- 
tions, 



(v\k) 



sin(fcr) 
kr ' 

the T-matrix equation (|25[) becomes, 

d 3 p v k pT pk , (E) 



T kk ,(E)=v kk , + 



(2tt) 3 E-p 2 /2M' 



(36) 



(37) 



The matrix elements of the delta shell potential are sep- 
arable: 



Vkk' 



sin(fcro) sin(/cVo) 
vo — ; n = v ou k u k * 



thus 

T kk XE)=v^ 



krn 



1 - v 



k'r 
3 d 3 p 



/,, (2tt) 3 E-p 2 /2M 
On shell (k = k' and E = k 2 /2M + in), 

v sin 2 (fcr ) 



(38) 



(39) 



T kk (k 2 /2M + in) 



k 2 r 2 + iMkv (1 - e 2ikr °) /4tt' 



(40) 



Comparing with Eq. (|27j) we find with some trigonometry 
that the s-wave phase shift for the delta-shell potential is 
given by 



2irkr 2 

cot<5(fc) = — — - -[1 + cot 2 (fcr )] - cot(fcr ). 

Mvq 



(41) 



Similarly in a sphere of radius R, with the s-wave func- 
tions 



1 sin(£; n r) 



V^kR 



(42) 
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as a basis, where k n = rnr/R, the T-matrix equation ([25 
becomes, 



1 



1=1 



E — E, 



Tin(E); (43) 



the matrix elements of the delta-shell potential are now 

vo 



2nRr 2 



sin(k m r Q )sm(k n r ). 



The solution of Eq. (|4U)l is jl4j 

OO . 2/7 \ 

w sm \hro) 



For E = fc 2 /2M, 



1 



E 



2-kRtI ^ E - kf/2M 



(44) 



(45) 



E 



sin 2 (kir ) R . 



i=i 
so that 



T- 1 = v~ 1 



— sm 2 (kr )[cot(kR) - cot(fcr )l, (46) 
2k 



Mv o . 2 
2nkr 2 



sin 2 (fcro)[cot(fci?) — cot(fcro)] 



(47) 



Comparing Eq. (g7]) with Eq. (|3TJ), we find that at the 
pole of T mn , at eigenenergy E = k 2 /2M, the expected 
result 



cot(kR) = — cot S(k); 



(48) 



this model calculation directly reproduces Eq. (j4j (with 
k = k' v ) within the T-matrix formalism for the delta shell 
potential, avoiding the subtle issue of the cutoff A in 
Eq. d3U). 



VI. PARTICLES IN TRAPS AND RYDBERG 
ATOMS 

So far we have considered containers in which the po- 
tential is zero outside the range of the scatterer. The 
results may be extended to particles in a potential, V(r), 
such as in an atomic trap. If the scatterer is at point 
r s , the Green function for the particle in the potential 
satisfies the equation 



1 

~2M 



(V 2 + k 2 ) + V(r) - V(r s ) 



G fe (r) = 5(r a ), (49) 



where k 2 = 2M[E - V(r s )}. Provided V(r) - V(r s ) is 
small compared with other energy scales in Eq. (J49j) at 
distances from the scatterer of the order of the range of 
the potential, the arguments of Sec. Ill still hold, with 
the origin shifted to r s . Explicitly, the condition is that 
\V(r)-V(r s )\ < 1/M6 2 for |r-r s | < b, in addition to the 
usual one kb < 1, where b is the range of the scattering 



potential (not V^(r)). Equation (Tl4l thus holds for this 
more general case, provided the <f> n are taken to be the 
eigenfunctions including the effects of the potential V(r). 

The general result (|14p can be applied to calculate the 
effective potential between a Rydberg atom and a small 
neutral atom. In the absence of the neutral atom, the 
eigenstates <^ n (r) for the outermost electron are hydro- 
genic wavefunctions for the Coulomb potential V c (r) ex- 
erted by the ion, where we take the nucleus to be at 
the origin. We place the neutral atom at position r s , of 
magnitude much larger than ao (the Bohr radius). The 
eigenenergies E of the outmost electron of the Rydberg 
atom arc functions of r s , and thus act as an effective 
potential between the neutral and Rydberg atoms. As a 
consequence of the electron's interaction with the neutral 
atom, the wave function of the electron in the Rydberg 
atom acquires a form near r = r s that corresponds to an 
s-wave phase shift 5(k). Thus from Eq. (TT4"|) , the effective 
potential E(r s ) is determined by 



4tt 



cot 5(h) = lim 



E 



i 



1 <MrR(r s ) 

2m e E n - E(r s ) 4?r|r - r s 



(50) 



where k = y/2m e [E — V c (r s )} is the scattering wave- 
length and m e is the electron mass. The sum over n 
includes both the discrete and continuum states. Solving 
this equation is an alternative route to the results of Ref. 
5] for the interaction of a neutral atom with a Rydberg 
atom. When applied to the problem of two particles in a 
harmonic potential interacting via a short-range poten- 
tial, Eq. (|50l) reproduces the result of Ref. [6j . 



VII. ENERGY SHIFTS IN MANY-PARTICLE 
SYSTEMS 

We give now an example that allows one to identify the 
quantity in the many-body problem that plays the role 
of the phase shift in two-body scattering, and that shows 
how the effects of the medium enter the scattering pro- 
cess. This section may be regarded as an illustration of 
the basic formalism for describing thermodynamic prop- 
erties of a many-body system in terms of fully renormal- 
ized single- and many-particle propagators given by De 
Dominicis and Martin [15| and reviewed in Ref. [16[ . We 
employ the delta shell model and calculate its effect in a 
dilute binary fermionic gas within the ladder approxima- 
tion [It], EE] ■ We assume that the system contains equal 
numbers of fermions in internal states |1) and \2) of equal 
mass m, with chemical potential /i. 

Differentiating the grand potential per unit volume, f2, 
with respect to the coupling constant vq gives 

on 

dv 



1 



d 3 p S(p - r ) (4 (p) (0)^2 (0)^i (p)) 



1 



B(p,u) 



pv 



v B(p,u))' 



(51) 
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where r\ is a positive infinitesimal and two-particle prop- 
agator is 



1 



P V Tt ( Z ~ £ p/2+q)( W - Z - e p/2- 



d 3 



(27T) 



9 2 

aV 



1 - 2/(e p/2+q ) 

2/i — p 2 /4m — q 2 /m' 



(52) 



e q = q 2 /2m — /j, to is a bosonic Matsubara frequency, and 
2 a fermionic one, / is the Fermi distribution function, 
(3 = 1/T with T the temperature (we set the Boltzmann 
constant to unity), and V is the volume. Here we ap- 
proximate the single particle propagators by those for 
free particles. In the low density limit, at zero tempera- 
ture, 1 — vqB(0, lu) is proportional to the inverse T-matrix 
for the separable potential and is defined in Eq. ([3T)|) . 

Integration of Eq. (|5Tj) with respect to vq gives the 
change in the grand potential due to the interaction: 



An 



-y 



d 3 p 



e w "ln [t( P ,u)/v ] 



where 



Vo 



1 - v B(p,uj) 



(53) 



(54) 



is related to the T-matrix in the medium by T^y {p, lu) — 
UkUk't{p,uj), with Uk = suikro/kro. 

We convert the sum in Eq. (|53"|) into an integral, choos- 
ing the branch cut of the logarithm to lie along the nega- 
tive real axis. We assume that 1 + mvo/4arro > so that 
there is no two-body bound state, and that the tempera- 
ture exceeds the BCS transition temperature. Then the 
only singularity of i(p, lu) is the branch cut along the real 
u> axis, so that 



d 3 p 
(2^)3 



dto 



f B (uj)5 m (p,uj) 



(55) 



where fs is the Bose distribution function. The quantity 



(56) 



6 m (p,u) = arg(t(p, u + if})) 



is the analog of 5 in the two body problem (cf. Eq. (|27[) ). 
Therefore 

An = J til) 3 (2tt) 3 ^( £ Pi)-f (^^(Pi'P 2 )- ( 57 ) 

By the theorem of small increments, small changes in 
thermodynamic potentials keeping the natural variables 
appropriate to the chosen potential are the same (l9j . 
Consequently, Eq. ([57)) gives the change in the energy 
density for fixed particle number and entropy per unit 
volume. Equation (1571) has the usual form for the change 
in the energy density for a two-body interaction, with an 
effective two-body interaction given by 



Weff(P!,P2) 



MlPi +P2|,e Pl +£ P2 ) 
IraBflpi + p 2 |, e pi + e P2 + iv) 



(58) 



In low density limit, — Im_B(|pi + pa|,e Pl + e P2 + irj)/ir 
reduces to m|pi — p 2 |/87r 2 , the density of states for a 
pair of atoms with total momentum |pi +pa|. For small 
|Pi — P2 1 , the phase shift is given in terms of the s-wave 
scattering length by S — > — |pi — P2|a s /2. The factor 
of two arises because in the problem of scattering from 
a rigid scattering center considered earlier, we worked in 
terms of k, the momentum of the scattered particle, while 
in the two-body problem, the corresponding quantity is 
|p x — p 2 1/2, the magnitude of the momentum of one of the 
particles relative to the center of mass of the pair. The 
effective interaction in momentum space then reduces to 
the well-known result 47ra s /m for two particles of equal 
mass (see Eq. (|17|) and the remarks following). 

The many-body calculation above shows that the total 
energy change is linear in a generalized phase shift that 
includes the effects of the medium on the scattering pro- 
cess. By contrast, Eq. (|T4"|) shows that in an arbitrary 
container, the energy changes of one-particle eigenstates 
in the presence of a short-range scatterer are generally 
not linear in the s-wave phase shift, and one similarly ex- 
pects the energy eigenstates of two or more particles with 
binary interactions in a container of arbitrary shape not 
to be linear in the phase shift for the two-body interac- 
tion. How does one account for this difference? The es- 
sential difference between the two situations may be seen 
by noting that the energy of a (non-interacting) Fermi gas 
in the presence of a central scatterer is a sum of contri- 
butions from a large number of individual single-particle 
energy levels. When the number of particles N is large 
compared with unity, the dominant contribution to the 
total energy is insensitive to the behavior of any partic- 
ular energy level. This is true even though the single- 
particle energies are not linear functions of 5, but exhibit 
the undulating behavior shown in Fig. [TJ Contributions 
beyond the leading term are, however, sensitive to the 
shape of the container and are responsible for, e.g., shell 
structure in nuclear binding energies. Similarly, in the 
many-body problem with two-body scattering described 
above, the total energy is insensitive to the boundary 
conditions a single-particle state satisfies on the walls of 
the container. 



VIII. CONCLUDING REMARKS 

One of the main results of this paper is a generaliza- 
tion of an expression due to Luscher [To| that expresses 
energy shifts of states in terms of scattering properties. 
In particular, we have derived a general expression, Eq. 
(jT4j) for the shift of the energy of the state of a single par- 
ticle in a container of any shape external potential when 
a scatterer of short range is added. In addition, a related 
result for a particle in a potential that is slowly varying 
in space was derived, Eq. (|50|) . The numerical results 
obtained for a particle in a cubic box show that the ener- 
gies of low-lying states are a rather complicated function 
of the phase shift. However, the functional form is con- 
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strained by the fact that when the phase shift is reduced 
from S to 6 — n, the state of the system is transformed 
into the next higher state for a phase shift 6. 

In Sec. VII, we considered the many-body problem and 
showed that, for a model system of two species of fermion 
interacting via a two-body potential, the shift in the en- 
ergy may be expressed as a linear function of an effective 
interaction proportional to a generalized phase shift that 
takes into account the effect of the medium on the scat- 
tering process. The reason for this simple behavior for 
the many-body problem is that the free energy is a ther- 
modynamic quantity, and this is sensitive mainly to the 



volume of the container, but insensitive to the details of 
the single particle wave functions, which depend on the 
shape of the container. 
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